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Abstract 

When Backlund transformations are applied repeatedly to any so- 
lution of the nonlinear equations describing electro diffusion through a 
liquid junction separating two infinite well-stirred layers, they give rise 
to quantized ionic fluxes across the junction. A particular exact solu- 
tion is shown to imply that these fluxes consist of quanta of electric 
charge. 
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1 Introduction 



Since the pioneering works of Nernst [T] and Planck [2], transport of charged 
ions across liquid junctions has played a fundamental role in a variety of 
natural systems. An extension of the Nernst-Planck model by one of us [31 E] 
incorporated the effect of the electric field that develops within a junction 
in response to diffusional separation of ions carrying different charges. The 
resulting nonlinearity is essential to what follows. 

In its simplest form, the model deals with steady one- dimensional trans- 
port of two ionic species with equal and opposite charges, across an infinite 
slab occupying < x < 5. The governing system of differential equations is 

c+ \ x ) = (ze/kT) E{x) c+(x) - $+/£>+ , 
c _ \ x ) = -(ze/kT) E(x) c_(x) - , 

E'(x) = (Anze/e) [c+(x) - c_(z)] (1) 

for < x < 5. Here c±(x) denote the concentrations of the two ionic species, 
$-t their steady (constant) fluxes in the +x-direction, D± their diffusion 
coefficients, and z their common valence, while E(x) denotes the induced 
electric field, k Boltzmann's constant, e the electronic charge, and T the 
ambient absolute temperature within the solution in the slab. An important 
auxiliary quantity is the electric current-density 

J = J + + J- , J± = ±ze $± . (2) 

Elimination of c + (x) and c_(x) from ([1]) results in a second-order nonlinear 
equation for E(x) that can be transformed into a Painleve equation of the 
second kind [S3 that has been studied extensively in this context [01 Ej. Here 
we prefer to keep to the set ([Tj and focus on the behaviour of all the physical 
variables. 

A feature [61 U\ E] of the system ([1]) is that it admits an auto-Backlund 
transformation [9 J B and its inverse B~ x . Thus, given any one solution 

S<°> = (cV 0) (x), cL 0) (z), E^(x), $i 0) , $L 0) ) , (3) 
a second solution is given by 

S (i) = g( 5 (o)) = ( c «( x ) ) c W( a .) > $ W ) (4) 
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where 

c+\x) = c_ (x) 



e ekT ( $f 

2nzeD + M( x ) + 2irz>e* Id +c <°>( 



x 



2&T $ 



(o) 



= 2$ (o) + ( jD+ / jD _)$(°) } = _p_/D + )$W , (5) 

and a third solution is given by 

= 5-1(5(0)) = (^(x), cL _1) (x), £ (_1) (a0, *\ , (6) 

where 



(-1)/ n = (o)/ e efcT / $L 0) N 2 

c ( : 1 )( x ) = c W(x), E(- 1 )(x) = -E(°)(x) 



ze£L c (0) (x) ' 

^l- 1 ) = 2$ (0) + , ^ = -(£>+/£>_)$<? . (7) 

Repeated application of B and £> _1 on any given seed solution generates 
a sequence of solutions 

S (n) , n = 0, ±1, ±2, . . . (8) 

that is doubly infinite in general. 

Considering repeated transformations as in (J5J), ([7]) and (|S]), we deduce 
that 

$W = ( n + + n( J D+/D_)$ (0) , 

$W = _( n - i)$(_ 0) - n(D_/D+)$ ( ® , (9) 



so that 



jf) = ( n + 1) Jf - n(D + /D_) J {0) , 
J W = -(n - 1) Ji 0) + n(D_/D + )jf , (10) 
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and hence 



jW = j(o) +nAJ) AJ = 5e(£> + + DJ) < + -jy- > , (11) 

where 

= Jf + Ji 0) = ($f - $i 0) ) . (12) 

This mathematical phenomenon has been called 'Backhand flux-quantization' 
[8]. It is associated with any solution of (JTJ taken as seed, in particular 
any exact solution, whether known [7] or as yet unknown. Physically the 
phenomenon may be related, for example, to discrete fluxes such as occur in 
electrochemical action potentials across excitable membranes [TU] in plants 
or animals. 



2 Quantal structure underlying solutions 

A notable feature of ( TTTj) is that increments of fluxes pertaining to successive 
values of n are all equal and all determined at n = 0, despite radical differ- 
ences between the solutions of ([[]) associated with those fluxes at different 
n values. What common element of such disparate solutions can engender 
consistency with equality of the flux increments? A special exact solution 
implies that the common element is the quantum of electric charge, as we 
now show. 

Consider the exact solution of ([[J with 

cf (x) = cL 0) (a) = c + (cx - c )x/5 , £ (0) (x) = , 

*® = D + (c -c 1 )/6, $L 0) = £>_(co-ci)/5. (13) 

Here Cq > c\ > are constants. This is a special case of a solution given 
by Planck [2] for the Nernst-Planck model, one that happens to provide 
also a solution of ([I]). It describes a steady-state situation where there the 
electric field vanishes throughout the slab, so that the ionic transport is 
purely diffusive, and charge neutrality is maintained throughout, including at 
the two faces. From this seed solution other exact solutions can be derived by 
Backhand transformations, as in (SD dHJ) • To simplify the discussion, suppose 
that the diffusion coefficients are equal: D + = D- = D. Then 

$ (0) = $ (0) = _ ^ 
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and hence (TTTT) and (fT2"j) give 

j(o) = _ J (o) = ^ ejD ( Co _ Cl )/ 5; 

J (0) = 0, AJ = 4^e J D(c -c 1 )/5. (15) 

It follows from (fT4"|) that in a time r, the number of positive ions (equal 
to the number of negative ions) flowing in the +x-direction (down the con- 
centration gradients) across an area A perpendicular to the x-axis is given 
by 

n+ = n_ = <S> { ± ] At = D(c - Cl ) At/ 8 . (16) 

Consider the situation where n+ = n_ « 1. Each of these two ions performs 
a random walk across the junction in a time 

t^5 2 /2D. (17) 

With this value for r, setting n + = n_ w 1 in ffT6l) means choosing A so that 

Aw2/(co- Cl )<J. (18) 

Then it follows from (IP21) and (|T5|) that this single positively (resp. nega- 
tively) charged ion carries a charge ze (resp. —ze) down the concentration 
gradients, and the nett charge across A is zero, consistent with the vanishing 
of the current density as in (fl5|) . Note that for these values of A and r, 

AJAr = 4ze. (19) 

Proceeding from to solutions 5 ( ' n \ n = ±1, ±2, we see from 
(fTTj) that AJ is unchanged despite changes in concentrations, appearance of 
electric fields, and loss of electroneutrality within the junction as evident in 
(jSD and (j7j). Thus electrical charges transferred in time r are, from (fTTj) 
and (USD, 

Q (n) = nAJAr = 4nze, n = 0, ±1, ±2, . . . , (20) 

so that the quanta implied here are quadruples of ze. [When D + ^ D_ the 
same incremental charges are transferred across the area A in the time 
t' = 2t + t_/ (r+ + r_), where t± is defined as in (IT7|) . with D-t replacing D.] 
To elucidate this result in terms of transfers of positive and negative 
ions, consider the Backlund-transformed solution as in (jSj), with Planck's 
solution f[T5|) as <S(°) . Keeping the same values for A and r we have from (fTTj) 

J (1) Ar = ( J (0) + A J) At = Aze . (21) 
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Furthermore, from ( llOp . 

J^At = 3ze , J W At = ze , (22) 

Thus describes a situation in which 3 positive ions are transported across 
A in the +x-direction in time r, while 1 negative ion is transported in the 
—x- direct ion. Note that transport in this situation is no longer by diffusion 
alone, as the electric field is now non- vanishing. 

The inverse-Backhand transformed solution <S( _1 ) describes the charge- 
conjugate situation, with 

J^At = -4ze, 



J { - 1] At = -3ze , J^At = -ze . (23) 

More generally, according to ffTUl) . ffTTl) and ffTIJl) . the solution describes 
a situation in which 

J {n) At = 4nze , 



jV*At = (2n + l)ze , At = (2n - l)ze , (24) 

for n = 0, ±1, ±2, . . . , corresponding to 2n+ 1 positively charged ions trans- 
ported across A in the +x-direction, and 2n — l negatively charged ions in the 
— x-direction. In this way, we see that Backlund flux-quantization directly 
reflects the quantization of electric charge on the two ionic species. 
Consider in more detail. From ( 1T3"|) and ([2]) we have 



c \Cq J 8 8nkTc 



where 



E ^ ( ze5 ) c + (ci - c )x/5 ^ 

with similar expressions for c[_ (x) and E^~^(x). The third (dimensionless) 
term on the RHS of (1251) typically has a value 1CT 10 or less [I], yet such 
terms give rise to the Backlund-generated exact solutions and S(~^ 
with Planck's solution (TT5T) as seed. While these derived solutions collapse 
onto the seed solution as eE^ 1 ^ (x) 2 /kTco goes to zero, the quantal structure 
(ED, (El, (El remains unaffected. 
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3 Historical note 

The theory of electrodiffusion was initiated by Planck [2] a decade before 
his quantization of black-body radiation [llj. The extensive developments 
of the latter were not slowed by mathematical complexities. In contrast, 
Planck's work on electrodiffusion had first to be supplemented by essential 
nonlinearities [3] H] before the pioneering ideas of Backlund [12] could be 
applied [6j [7J [8] to reveal at least a fraction of the wealth of exact solutions 
of (pQ), based only on the particular seed solution (Tl3|) . but sufficient to imply 
the quanta of electric charge, as in (1201) and (1241) . Despite the passage of 
time, these surprising products of the nonlinearity of the continuum theory 
(PQ) may be regarded as Planck's other quanta. 
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